Trapped slender vortex filaments 1 
The statistical mechanics of nearly parallel vortex filaments confined in the unbounded 

plane by angular momentum, first studied by Lions and Majda (2000), is investigated 
using a mean-field approximation to interaction and a spherical constraint to develop 
an explicit formula for the mean square vortex position or length scale of the system, 
R, verified with Path Integral Monte Carlo simulations. We confirm that 3D filaments 
resist confinement in a different way than 2D point vortices and that this results in a 
profound shift at high-densities for the length scale of quasi-2D versus strictly-2D models 
of vorticity fields in which angular momentum is conserved. Our analytical results corre- 
spond well with those of the Monte Carlo simulations and show a 3D effects contributing 
significantly to determination of the length scale. 
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1. Introduction 

Statistical mechanics is a way of calculating the macroscopic properties of matter from 
the probabilistic behavior of its microscopic components. Rather than solving the Navier- 
Stokes equations explicitly in time, a statistical equilibrium approach for fluid flows aims 
to describe observable quantities by averagin g over "microstate s " or s tates that account 



for a system's exact position in phase space (jMaida and WansI (j2006[ )). 

In fluid turbulence a statistical treatment is often preferable to obtaining direct solu- 
tions to the Navier-Stokes equations because of the inherent chaos/complexity of turbu- 
lent flows that makes it impossible to model every trajectory at high Reynolds numbers. 
This nonequilibrium statistical description of turbulence can be replaced by a statistical 
equilibrium approach under specific conditions such as nearly inviscid quasi-2D flows 
where a separation of time scales is valid. Familiar arguments in favor of a separation 
in the time scales of energy and angular momentum transfer at flow boundaries versus 

viscous dissipation are based on the tea cup paradigm. The "universal equilibrium as- 
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sumption" , that the time scale at which turbulent features form is much small er than 



the tim e scale for viscous decay, justifies studying turbulent flows in equilibrium (jChorin 



( 19941 )). Although this justification can be considered tenuous and subject to mathemat- 
ical limits, experimental observations provide a validation for the statistical approach. 

Most work on fiuid flows in statistical equilibrium is on how large scale structures 
appear in vortic ity fields of 2D rotating ideal fluids such as Onsager's Point Vortex Gas 



(jOnsageii (jl949l )). Many significant results have come out of this approach, but research 
has confirmed few 2D results for nearly-2D or quasi-2D models. In many cases, 2D models 
approximate a quasi-2D reality. This paper is concerned with when this approximation 
fails and how the statistics of quasi-2D models depart from those of fully-2D when it 
does fail. 

Nearly parallel vortex filaments (Figure [2]) are one of the simplest model for discrete, 
C|uasi-2D vor t icity. They appear in an incompressible, nearly inviscid Navier- Stokes flow 



( Klein et al 



(1995!)) as well as other physical systems such as in rotating s uperfluids 



and Bose-Einstein conden s ates. F rom theories for single vort ex fllaments by 



(1972) 



Callegari and Tind (|l978l ). 



Klein et al 



ments. 



Ting and KleinI (|l991l ). and 



Hasimoto 



Klein and Maidal (|199lD . 



19951) develop ed the first rigorous model for interacting nearly parallel fila- 



Julien et al 



(|1996f ) have observed these nearly parallel vortex filaments in astro- 
physical simuations. 

Our statistical de s criptio n of nearly parallel vortex filaments derives from the work 



of iLions and Maidal (j2000[ ). In their paper, they developed a description of the time- 
dependent statistics, namely the density, p{t,r), where r e 3?^, of the filaments for 
N —f oo, given tha t energy and angular momentum are conserved, using the model 



of 



Klein et al 



( 19951 ). Filaments are also periodic in the z-direction with period L. 



They introduce a Gibbs probability measure for the vortex filaments in the unbounded 
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Figure 1. Shown here in top-down projection, nearly parallel vortex filaments, at low density 
and high strength of interaction, are well-order ed into a 2 - D tri angular lattice known as the 



Abrikosov lattice from type-II superconductors IjAbrikosovl ([1957|)). This figure shows how the 
quasi-2D model is essentially a 2-D model for these parameters. 

plane, confined by angular momentum. This distribution models a system that is not com- 
pletely isolated, a reasonable assumption in geophysical and other macroscopic flows, in 
that a quiescent fluid surrounding the filaments acts as an energy and angular momentum 
bath. The system is allowed to exchange angular momentum and energy with the envi- 
ronment and so angular momentum and energy fluctuate. A justification, a posteriori, of 
the Gibbs' probability measure that is canonical, rather than microcanonical, in energy 
and angular momentum comes from results indicating that vortices are confined to a 
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Figure 2. Nearly parallel vortex filaments are nearly parallel to the z-axis in an asymptotic 
sense that their deviation from straight is a small parameter. They are infinite in length but 
have a period L. Here there are 50 filaments. 

compact domain in the plane, and the heat bath exists outside this domain. In formal 
notation, a canonical Gibbs' measure means that if a state s has energy Es and angular 
momentum Is, then the probability of s, 



P{Es,L)o,e-^^^-'^'% 



(1.1) 



where /3 and fj, are constants. (For historical reasons (3 is known as inverse "temper- 
ature" and fi as "chemical potential" although they are not directly related to their 
molecular equivalents because the "molecules" here are vortices . ) Sinc e J2s-P(^s,Is) = 



l,P{Es,Is) = Z-^e-^^^^-^'^% where Z = Y^^e-l^^^-^'^' jchorinl (|l994h ). Lions' and Ma- 



jda's rigorous treatment of this statistical ensemble leads to a non-linear Schroedinger 
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equation governing p. Their approach has been criticis ed as being physically inconsistent 



because it admits states that are not nearly parallel (|Berdichevskvl (|2002l) ). This criti- 
cism is valid. However, a careful choice of parameters /3, p,, etc. causes these disallowed 
states to have such large energies that their probabilities are negligible. For such param- 
eter ranges, we can assume that the Lions-Majda model is physical. Our Monte Carlo 
simulation naturally rejects states with non-nearly-parallel lines without additions to the 
model. 

The conservation of angular momentum is key in this study because it in troduces a nat- 



ural le ngth scale to the confined system of vortices in the unbounded plane (IDiBattista and Maida 



20061) 



( 2001 ) ). It also plays an important role in statistical mechanics on the sphere 



LimI (|2006l) 



Maida and Wane 



Lim and NebusI (|2006r )l. Periodic boundary conditions are more com- 
mon, but these enforce a length scale (the period) artificially on the plane that can affect 
the statistics. Although a pre-determined length scale is reasonable in some cases, for 
many others, it is not. Moreover, the statistics for periodically bounded ensembles is 
quite different from unbounded systems that are confined by angular momentum or a 
harmonic trap as in the case of Bose-Einstein condensates. As the reader will see in this 
paper, our results, which focus on determining the length scale both analytically and 
computationally, depend upon the angular momentum confinement. 

Our goal is not merely to demonstrate that nearly parallel vortex filaments behave 
differently from 2D point vortices. Nearly parallel vortex filaments at low-density /high- 
straightness (large /3 or low-temperature) behave exactly like 2-D point vortices (Fig- 
ure [T]) . That their behavior at high densities and levels of curvature (small (3 or high- 
temperature) ought to be different than that of 2D point vortices is obvious. On the 



other hand, despite the rigorous work of Lions and Majda, the qualitative statistical 
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differences between nearly parallel vortex filaments and strictly-parallel filaments or 2D 

point vortices are unknown. 

Our hypothesis is that at high-densities point vortices and nearly parallel vortex fila- 
ments have qualitatively diff erent behavior. The len gth scale of point vortices collapses 



with increasing temperature (jLim and AssadI (|2005r )). We hypothesize that filaments re- 



verse this collapse (Sec.|6]). Similarily to how stars in globular clusters resist gravitational 
collapse through motion, we posit that filaments resist collapse through curvature (which 
in top-down projection appears like Brownian motion) and that the outcome is a reversal 
of collapse as temperature increases. This idea has profound implications for the natural 
length scale of 2D versus quasi-2D models. 

We develop a simpler approach than Lions and Majda — one that is motivated and 
justified by their rigorous mean-field result — that we hope generates a more intuitive 
understanding of what happens when the density and curvature of the filaments increases 
and interaction becomes more three dimensional. Our method is two-fold: 

(a) We develop an explicit formula for the second moment, R^, (i.e. the square length 
scale) of the equilibrium vortex density distribution, p(t — > cxd, r), using a simpler mean- 
field approximation to interaction than Lions-Majda. Our mean-field approach assumes 
that the interaction due to iV — f filaments on one filament is similar to how a perfectly 
straight filament fixed at the origin with strength — 1 would affect the the center-of- 
mass of that filament. This assumption requires that the filaments be fairly uniformly 
distributed, whereas Lions-Majda do not require uniformity, ft is a "particle interacting 
with a center-of-mass" rather than "particle interacting with a density field" assumption. 
We then constrain the filament's planar position such that, 







rfr|1'(T)|" = Li?^ (1.2) 
where vE'(r) — x{t) + iyir) represents the curve of a filament in complex notation. 
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This constraint enters the Gibbs distribution (Eqn. II. ip as an exact (micro-canonical) 
conservation law. These two assumptions, discussed in detail below, allow us to derive 
the formula: 

2 _ p'^aN + \/;a4^^2]v2:jr32^ 



(1.3) 



where a is related to core structure. (Larger a means straighter filame nts.) (Sec t ion [4 ) 



(&) To confirm our formula, we perform Path Integral Monte Carlo ijCeperlevI (|l9951 )) 
on the original statistical system using piecewise linear approximations to the filaments. 
These simulations also allow us to de t ermin e that our results are physical, addressing 
the point of criticism bv iBerdichevskvi (|2002f) . (Section [S]) 



Our statistical derivation not only answers the questions above but gives an extremely 
close prediction of length scale, R, and allows us to explain the qualitative differences 
between a strictly-2D and a nearly-2D model, giving new insight into the role of 3D 
effects in the onset of turbulence in quasi-2D models. 

The paper is laid out as follows: we discuss the problem i n greater detail, i nclud- 



Lions and Maidal (|2000l) 



ing the 2D Onsager Gas (SJ2]), go over the model presented in 
(® 1 present our mean-field derivation for R^ O, a description of Path Integral Monte 
Carlo (il]), and a comparison of Monte Carlo results with the mean-field formula and a 
verification of the assumptions of the model (^H). We discuss and conclude (33 



2. Problem 

The low-temperature (asymptotic as /3 — s- oo) statistics of strictly-2D point vortices 
in the unbounded plane with angular momentum conserved are well understood. The 
Point- Vortex Gas Model has the Hamiltonian, 

HlF = -Y,X^,log\^.~^,\^ (2.1) 

j>i 
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where Xi and Xj are the circulation constants for point vortices i and j and and 

are their positions in the complex plane. This Hamiltonian derives from the 2D Euler 

equations in vorticity form given that w, the field, has the form 



where is a position in the complex plane. (jChorinI (|1994[ )). The probability distribution, 

exp {-PHn - (J-In) 



P. 



N 



where 



and 



Zn = j d^i--- J d-^N exp {-f3HN - ^iIn) 



(2.2) 



(2.3) 



N 



/jv=^A,|vI/,|2, (2.4) 

i 

represents a non-isolated system in which energy and angular momentum are exchanged 
with the environment at some level of fluctuation determined by the "inverse tempera- 
ture" , (3, and the "chemical potential" , fi. Note that it is possible for the temperature 
here to be negati ve but not the chemica l potential for angular momentum. 
In their paper. 



Lim and AssadI (|2005f ) show variationally that the mean square vortex 



position (variance) for 2D point vortices, 
where the average is over P/v, has a nice formula 



R 



2D 



(2.5) 



(2.6) 



where fl = A,; is the total circulation. In our case, we define A^ = IVi, so Q, ~ N and 
i?2£) = NP/{Aii). Moreover, they were able to show in Monte Carlo simulations of an 
ensemble of 1000 point vortices that the distribution of vortices is almost uniform and 
axisymmetric, meaning that the probability distribution of vortices is nearly a perfect 
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cylinder, suggesting that R^j^ is the only value we need to know to determine the entire 
distribution. 

These results provide a tantalizing starting point for our investigation of nearly parallel 
vortex filaments where, obviously, at some level of density and fluctuation there exists 
different statistical behavior. Moreover, the discovery of a uniform distribution at low- 
temperature, we hope, carries over to nearly parallel vortex filaments. 

The low-temperature formula for i?2D ' Equation 12.61 decreases with /3, meaning that 
the radius decreases as the "temperature" , related to the speed of motion of the vortices, 
increases. At some point, of course, /3 is so small that the distribution becomes solely 
a function of angular momentum, clearly a normal distribution with variance, R^jj = 
l/(2/i). The radius does not decrease to a point, but it never increases, certainly not in 
the low-temperature regime. 

The question we address in this paper concerns the statistical distribution of nearly 
parallel vortex filaments at low to moderate temperature, not high-temperature. There- 
fore, our derivations in Section |4] assume that (3 is large enough that the fluctuations 
in interaction energy and angular momentum do not play a major role in the statistics, 
only their mean- values. The role of the Monte Carlo comparison in Section[6]is to justify 
this assumption. Specifically, we are concerned with the mean square vortex position, 
i?^, which we define below in Section [4] 

We answer the questions: (1) do nearly parallel vortex filaments resist and reverse 
planar collapse as temperature increases (an important question if one would like to 
contain a column of rotating fiuid) and (2) can we see this behavior without violating 
assumptions of nearly parallel, i.e. do the filaments have to be too curvy before they 
begin to reverse collapse? Neither answer is obvious. 
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3. Model 



Our model properly derives from the paper of 



distribution that they introduce. Although iKlein et al 



Lions and Maidal |2000f) and the Gibbs 



19951 ) derived the equations of 



motion, we make use of the statistical framework presented in the later paper and also 
we rely heavily on their broken-segment model, not only in the Monte Carlo simulation 
where it is necessary to discretize the z-axis but also in our mean-field approach. 
We define nearly parallel vortex filaments. 

Definition 3.1. Nearly parallel vortex filaments are smooth curves that we can repre- 
sent with a complex parameterization ^,;(t) where — Xi{T) ~\- iyi{T) and r S [0, L). 
They have a special asymptotic form, Given that i G 0(1), if we take any two values of 
T, Tq and Ti such that Tq < Ti, and let At = Ti — Tq such that At G 0{e) where e <C 1, 
then for any filament i, the amplitude |5'i(ri) — ^i(to)| € O(e^). 

In words this means that, for a small rise of length e in the filament, the amplitude 
must be on the order of e^. This assumption guarantees a certain degree of straightness 
in the filament that allows for the derivation of the quasi-2D equations of motion. The 
other asymptotic assumption is of the vortex core-size, h, which has the property /i <C e. 
In this model we assume that the vortices have no cross-section, i.e. the vorticity field, 
uj, has the form, 

u; = ^<5(vI/(r)-^P,(T)), (3.1) 

i 

where ^'(t) = x{t) + iy{T) corresponds to the 3D Cartesian position (x(t), ?/(t), t). 
Equilibrium statistical mechanics is trad itionally concerned with conserved quantities 



in a Hamiltonian system. The Hamiltonian 



Klein et al 



( 19951 ) derived for nearly parallel 



vortex filaments has the form. 



Hn = OL 



•'^ i=l 



5*.(t) 



dr 



N N 

drJ2Y.^og\^^iT)~^Jir)\, (3.2) 

2—1 J>2 
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where a is the core-structure constant in units of energy/length. This Hamiltonian re- 
sembles the 2D point vortex Hamihonian in that the interaction is logarithmic in the 
plane. For two filaments i and j, only points in the same plane, the same value of r, 
interact. Two points at different values of r only interact if they are both on the same 
filament. The first term is a local self-induction term that essentially causes the filaments 
to wriggle in top-down projection like a particle under Brownian motion. The additional 
conserved quantity, angular momentum, is given by, 

'0 



In^Y. rdr\^dr)\'- (3.3) 

, Jo 



Lions and Maidal ([20001) introduce a broken segment model in Section 2.2, Equation 
2.20 of their paper. For each i, 'l'i(r) is piecewise linear with AI segments. Each vertex 
or "bead" is at a multiple of (5 = L/M . Therefore, we define 7/'i(fc) = — 1)<5), 

and the curve for filament i is given by the vector — (■(/'i(l), . . . ,ipi{M)). With this 
representation we rewrite the conserved quantities, 

H^iM) = a E E iMl(^±lp:ii^ - E E E ^ l^^(^) - (^)l' (3-4) 

i=l k=l i=l j>i k=l 

N M 

J^(M)=^^<5|V.(fc)r (3.5) 
i=i k=i 

The Gibbs probability measure for a state s = (\l/i, . . . , ^jv) is then 

Pis) = expi-PH^iM) - iilNiM)), (3.6) 
where the partition function has the form, 

Zn{M)= I d^i- ■ ■d^Nexp{-f3HN{M) - filN{M)). (3.7) 
We note that while it is possible to simulate the system with Mo nte Carlo it is not possible 



Lions and Maidal (|2000l) 



to solve explicitly for Zn for any value of AI. In their paper, 
go on to derive a non-linear Schroedinger equation that can give approximate values for 
P(s). Their PDE captures a great deal of the statistics, but they do not provide any 
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explicit formula for the length scale R, which is our goal. We refer the reader to their 

paper to learn more about their derivation and the model. 

4. Free Energy Theory 

4.1. Free Energy of Most- Probable Macrostate 

Given a functional for the free energy for a system, it is possible to solve for the 
statistics of the most-probable macrostate by minimizing F with respect to the desired 
statistic. In our case the statistic is defined as: 

Definition 4.1. The mean square vortex position, 

R^^{L-'N-'J2 / ^^I^^WI')' (4-1) 
where ^^(t) = Xi{T)+iyi{T) is a complex number representing the 2D position of filament 
i at z — T and the average (*) is with respect to the 3D Gibbs probability measure in "^i, 
Pn — ^ d^i ■ ■ ■ d^ isi e'Kjp{-' j3Hfq — ^jlIn), where Hfq is defined by Equation \3.S\ In 
by Eauation W^ 

The difficulty lies in deriving that functional. In systems such as ours, with fixed volume 
iy — 3?'^) the Helmholtz description of free energy is appropriate, F = U — TS, where 
U is average energy, T is temperature, 5* is entropy. 

In our system we write the free energy as follows: 

Fn = {Hn) + ^{In)~^Sn, (4.2) 
P P 

where S^Sn,U ^ (Hn) + ^{In), and T = i. 

While many statistical mechanics approaches are devoted to finding a functional for 
Sm, since and are usually known, there is another way to develop the free energy 
functional from the partition function Z^. The function Z^ can be considered a sum 
over microstates, Si, or it can be an average over a set of macrostates, Cj. Equation 12.31 
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is a sum over microstates. Alternatively, the partition function is given by the formula 

Zm = ^exp(-/3i/jvK-] - iilN[<yj])P[<yj], (4.3) 

j 

where Hj^i is the energy, Iiq is the angular momentum, and P is the probability for 
macrostate <Tj. 

Since SM[<yj\ = logP[crj], using the formula 14.21 we can say, 

3 

= ^exp(-/?FjvK-]). (4.4) 

3 

Because of conservation laws, we assume in physics that the most-probable macrostate 
or energy-state, j = m, has a probability so much larger than the probabilites of all other 
macrostates that sum contributions from other macrostates can be neglected and 

Zn = eM-PFN[a„^]). (4.5) 

Therefore, we have an equation for the free energy, 

Fat = -ilog^AT. (4.6) 

Even though this free energy is only the free energy of the most-probable macrostate, it 
can be considered the system's free energy. 

4.2. Mean Field Assumption 

As mentioned in Section [3l we cannot solve for Zjsi. We need to approximate it to derive 
the free energy functional. One way to do this is with a mean-field assumption. Because 
the logarithmic interaction in the Hamiltonian prevents us from evaluating the integral 
over microstates, Equation l2.3( we need to remove the interdependence of filaments upon 
one-another so that we can separate the integ rals over different filaments. In their paper 



this is the approach of 



Lions and Maidal ([20001), in that they develop a model of a filament 
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that interacts with a probabihty distribution, or field, rather than other filaments. Our 

approach is to take a cue from the physics of Newtonian gravity and assume that the 

filament interacts with the other filaments in the same way its center of mass would 

interact with an imaginary, "center of mass" fil ament. We mot i vate t his assumption 

from the Monte Carlo results on point vortices of iLim and AssadI (|2005h that showed a 



"fiat-top" , cylindrical probability distribution for vortices in the plane. A nearly uniform, 
symmetric distribution is crucial to the center-of-mass assumption. 

Given a point on filament i, rpiir), if the filament's center-of-mass interacts with the 
center of mass of the other filaments, the interaction potential, V^, simplifies: 

1 



2 



-(iV- l)-log|V',(T)| = -iV/41og|^,(T)|2 



'^dTmr)\^ 



= -iV/41og {^^^ 
= -^/41og^, 

where we can say for large N that iV — 1 ^ A^. This result makes the interaction a 
function of the angular momentum. 

The center-of-mass assumption liberates us in the evaluation of the partition function, 
because we need to consider neither the configurations of other filaments nor their density 
distribution, whereas Lions and Majda do take the latter into account. All we need 
consider is the angular momentum. 

4.3. Spherical Constraint 

Given the mean-field, center-of-mass assumption, we have a new Hamiltonian system 
governing the behavior of N independent, Af-segment filaments. The Hamiltonian reads 



16 
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i=l k=l 2 = 1 fe = l 



(4.7) 



2 5 ' ° L 

2 = 1 fe = l 

where H'fP{M) only differs from Hm{M) (Equation 13. 5p in the second term, and the 
new partition function reads 

Z]v"(Af) = y • • • d*jv exp(-/3i/™'(Af) - fiPrTiM)) 

dM'iexp(-/3i7r(M)-/i/r(A^))| , (4.8) 

because all the integrals are independent and equal. 

Returning to our derivation for the length scale, R, we take the simplest course possible 
to give a reasonable formula and impose a constraint on the filament, (dropping the 1 
subscript): 



/ dT\^{T)\'=R'=Y,M~'\m\ 



(4.9) 



which imposes an average distance from the origin over the length of the filament. We 
include this constraint as a micro-canonical or exact constraint into the partition function, 

N 

. (4.10) 



fjcm^sc 



This constraint eliminates fluctuations in interaction energy and angular momentum, 
leaving only fluctuations in local self-induction energy. This eli mination provides the 



Lim and AssadI ( 2005 ) 



best opportunity for comparison with the 2D length scale work of 
that also eliminates these two types of fluctuation. 

Eauation l4.10l is a multi-dimensional Gaussian integral with a spherical constraint — so 
called because it forces the vector to stay on the M — 1-sphere of radius \fMR. To 
evaluate it we turn to the spherical model of 



Berlin and Kac 



1952). 



The spherical model comprises a number of steps for evaluating integrals of the form 
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of Equation I4.10| beginning by putting the Dirac-delta function into integral form and 

ending with a steepest descent evaluation of the integral over ^. 

In integral (Fourier) form the Dirac-delta reads: 



-oo 27r ^ 



M 



MR^-Y,\m\' 



k=l 



(4.11) 



and allows us to combine the function in the exponent in Equation l4.10l with the exponent 
of the spherical constraint. (Whereas before a was used as a symbol for macrostate, here 
it is an integration variable.) 

4.4. Free Energy Derivation for R 

Now we can determine the free energy functional, which ought to be minimal under the 
constraints of the system, via steepest descent. Determining the free energy as a fmiction 
of the particular statistic, R, is more useful than having an equation for the partition 
function itself because we can determine R by minimizing the free energy w.r.t. it and 
then solving for R. 

Theorem 4.2. Given the partition function defined in Equation \4-10\ and the free 
energy in Equation\4-.2\ as M — s- oo, the value for R^ giving minimal free energy is, 



= —■ (4-12) 



Proof. The steepest descent method is an excellent way to obtain values for the 
minimum free energy. In employing steepest descent, we need to start with an integral 
of the form, J dxe^^^^^^^ . Then if -F[a;o] < i^[a::]Va: ^ xq, i.e. its minimum value is at xo, 

Jim M-Hog(^J dxe^^^^l^l^ ^ F[xo]. 

This works because, as M increases, the distribution that e^^-'^^I^l represents becomes 
narrower and focuses on xq until the distribution has zero value at all other x. 



18 T. D. Andersen and C. C. Lim 

To take this approach, we put the partition function in the appropriate form: If 

then 



2 



^ iiLM-^R^ - iG{R^ - y^M-\'il,{k)'^ , (4.14) 

is the non-dimensional free energy, where we have already applied the Dirac delta function 
to the interaction energy and the angular momentum. 

We can pull the constant terms out of the integral over ^. Because F is positive 
definite, under Fubini's theorem we may switch the integrals to obtain 

Z'^'^^M) = |e^/3V4iogK^-/.Lfl= ^ / '^^e-^'I^^lj'^ , (4.15) 



where 

fe=l \fe=l / 

is the free energy still dependent on "i. 

The interior integral needs evaluation. Let s = ia and 

/ioo T /> 
^ / d*e-^'W, (4.17) 

be that interior. Let us put F' in matrix form: 

F'[s] = -sMi?2+i^^t^^ + s^t^^ (418) 

where K = a(3/5, and the M x M matrix A has the form 

= = Ai^M = Am,1 ~ ~2' 

Aij = other 
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The integral in Equation l4.17l is Gaussian. We can evaluate it, knowing the eigenvalues 



of the matrix A. These eigenvalues have the for m Aj = 1 — cos(27 r (i — 1 



to the previous use of Xi as strength of vorticity) (|Berlin and Kad (|l952l ) 



) /M), (not related 



Lions and Mai da 



(|200a) ) and so 

Z'{M) = r |le^^«V*^n— W- (4.19) 
J-joo 27r Y * + 

We need to put Z'{M) back into the correct form for steepest descent. Following the 
example of Berlin and Kac, let s — K{rj — 1) then 

/ioo J 
^Ktt^Hv - l)-ie-*^i°gWe^'^/M, (4.20) 

where 

M 

= KR^{ri - 1) - M-^ log(?7 - cos(27r(i - 1)/M)), (4.21) 

i=2 

and rj > 1. 

We leave the i = 1 term out of the sum in f[r]] so that we can evaluate / further by 



taking the limit on the second term, 
m 

M- 



lim y^log(77 - cos(27r(i - 1)/M)) = — / duj\og{r) - cos{u))) 



i=2 

which gives 

M = KR^ri -1)-1-J^ diulogiv - cos{u)) 

= i^i?2(77-l)-log(77+(r;2-l)^). (4.22) 

To apply steepest descent, we determine the saddle point rj = rjo where f[ri\ has its 
minimum value, /[t^o]- Taking the derivative and setting it equal to zero, 

?f = KR''- ^ ^ = 0, (4.23) 

implies 



20 T. D. Andersen and C. C. Lim 

Having evaluated /, we can give an equation for the original free energy, 

F[7?o] = -NpL/AlogR^ + iiLR" + Mlog{K) - Mf[7jo], (4.25) 

and evaluate it as M — > oo. 

The term in the limit M log{K) , docs not depend on , and it is an unnecessary 
component representing the entropy of the broken filaments in the non-interacting case, 
and we drop it. Now we fill in the expressions for rio and K as defined above: 

lim M/[%]= lim KR^rjo - 1) - M log{vo + {vl - '^)^) 

M^oo M^oo 



hm MKR^ i , —Lr-- + 1 - 1 I -Mlogi J .„L.. + 1 ^ 



M^oo \ Y {KR^'f / I Y (iCi?2)2 KR^ 



= lim M^^R^ ij^^ + 1- 1 I -Mlog ( J— + ^ + ^ 

(4.26) 

which, because it is an energy for a filament, ought to be finite. The first term is the 
energy of the filament, Efii, and the second, the entropy, Sfu, and each by itself is finite, 
so we take each limit separately. 

The energy limit is simple to calculate using L'Hopital's rule: 



Efu= lim M^R^J 

M^oo L I Y {oipMR^f 



lim Am flo +1 



5^0 2S \\{apR^f J {apR^f 
L 

2a/3E2 



(4.27) 
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The entropy limit is equally simple: 



Sfii — lim Af loE 



L 

= lim — log 

.5^0 5 



lim L 



1 



1 



1 



s 



(52 



-1 



(a/3i?2)2 



a/?i?2 



;a/3i?2)2 



a/3i?2 

These two results imply that 



and 



lim M/frjol = ^ 



We minimize with respect to , 

dF 



and solve for i?2 



5i?2 



i?2 = 



= L/i — 



4i?2 " 2a/3i?4 



2a/3i?2 ■ 



= 0, 



iV/3/4±^(iV/3/4)2 + 4^^ 



2^ 



where we take the "plus" solution as giving physical results. 
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(a/3i?2)2 a/3i?2 

(4.28) 



(4.29) 



(4.30) 



(4.31) 
□ 



We note that the relationship between ^[770] and our previous Fjv in Equation 14.21 is 
PVio] — PFn, meaning that F[r]Q] is non-dimensional while Fn has units of energy. 
The resulting expressi on for the squar e lengt h scale is useful for comparison with 



the length scale result of iLim and AssadI (|2005f l because, if we take the limit 



lim 

CK — >00 



8a/?/i 4/i ' 



(4.32) 
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we get back the 2D poin t vorte x result for the length scale, which shows that our formula 



and the iLim and AssadI (j2005[ ) formula agree for perfectly straight filaments. 

However, for finite a the two formulae show a signficant difference. The 2D formula is 
linear in /3, our formula is non-linear. In fact, for decreasing /3, the sign of the slope of 



our formula changes at /3 = /3o, where 



f3i 



4m 



(4.33) 



That the system collapses and then starts to expand as "temperature", 1/(3, increases 
indicates a significant departure from the strictly-2D where the system size only collapses. 
In SectionlU we show that that Monte Carlo confirms this result and that the straightness 
assumptions of the model hold through much of the expansion phase. 



5. Monte Carlo 



Path Integral Monte Carlo methods emerged from th e pa t h inte gral formulation in- 



vented by Dirac that Richard Feynman later expanded (jZed (|2003f) ). in which particles 
are conceived to follow all paths through space. One of Feynman's great contributions 
to the quantum many-body probl em was the mapping of path i ntegrals onto a classical 



system of interacting "polymers" (jFevnman and Wheeled (|l94a )). D. M. Ceperley used 



Feynman's convenient piecewise linear formulation to develop his PIMC method which 
he successfully applied to He-4, generatin g the well - know n lambda transition for the first 



time in a microscopic particle simulation (jCeperlevI (|l995f l). Because it describes a system 



of interacting polymers, PIMC applies to classical systems that have a "polymer" -type 
description like nearly parallel vortex filaments. 

PIMC has several advantages. It is a continuum Monte Carlo algorithm, relying on 
no spatial lattice. Only time (length in the z-direction in the case of vortex filaments) is 
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discretized, and the algorithm makes no assumptions about types of phase transitions or 

trial wavefunctions. 

The Monte Carlo simulation begins with a random distribution of filament end-points 
in a square of side 10, and there are two possible moves that the algorithm chooses at 
random: 

(a) Moves a filament's end-points, and + 1). The index i is chosen at 

random, and the filament Vs end-points moved a uniform random distance. 

(6) Keeps end-points stationary and, following the bisection method of Ceperley (Fig- 
ure 131) , grows a n ew internal configuration for a randomly chosen filament 



(ICeperlevI (11995D). 

In each case, the energy of the new state, s', is calculated and retained with probability 

A{s -> s') = min{l,cxp (-[3[H%{M) - H%{M)] ~ fiir^iM) - /^(M)])} , (5.1) 

where s is the previous state. This effectively samples states from the Gibbs probability 
distribution in Equation l3.6l 

Our stopping criteria is graphical in that we ensure that the cumulative arithmetic 
mean of the energy, 

k 

ELm = ^^(^^) + f ^^'(^■')' (5.2) 

i=l ^ 

where Si refers to the state resulting from the ith move and k is the current move index, 
settles to a constant. The energy is almost guaranteed to settle in the case of the Gibbs' 
measure because of the tendency for the system to select a particular energystate (mean 
energy) and remain close to that state. Typically, we run for 1 million moves (accepted 
plus rejected) or 50,000 sweeps for 20 vortices. Afterwards, we collect data from about 
200,000 moves (1000 sweeps) to generate statistical information. 
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12 3 infinity 



Figure 3. The bisection algorithm works by bisecting the filament to sample point positions. 
First the center point is selected, then the two points half-way from the center to the end-points, 
then four more points, eight, and so on until some maximum number of points are sampled. 
These are snapshots of one filament at different steps in the sampling process. 

6. Results 

6.1. Comparison 

We simulated a collection of iV = 20 vortices each with a piecewise linear representation 
with M = 1024 segments and ran the system to equilibration, determined by the settling 
of the mean and variance of the total energy. We ran the system for 20 logarithmically 
spaced values of /3 between 0.001 and 1 plus two points, 10 and 100. We set a — 10^ 
(enforcing straightness) , fj. = 2000, and L = 10. Decreasing /3 simulates an increase in 
temperature, 1//3. We calculate two arithmetic averages: the mean square vortex position, 

N M 

i?i,c = (M7V)-i^^|Vi(fc)p, (6.1) 

i=i k=i 

and the mean square amplitude per segment, 

N M 

a' = (MiV)-i ^ Yl W - ^^(^ + 1)1'' (6-2) 

i=l fe=l 

where ipi{M + 1)= tpi{l). 




10-6 I I 
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Figure 4. The mean square vortex position, defined in Eauation l6.1l compared with Equations 
ll.3l and[ 2.6l shows how 3-D effects come into play around /3 — 0.16. That the 2D formula continues 
to decrease while the Monte Carlo and the quasi-2D formula curve upwards with decreasing /3 
suggests that the internal variations of the vortex lines have a significant effect on the probability 
distribution of vortices. 

Measures of the Monte Carlo Rf jc Equation 16. H correspond well to the 3D R^, 

Equation ll.31 in Figure |4] whereas, the strictly-2D i?!^), Eauation l2.6i continues to decline 

when the others curve up with decreasing f3 values, suggesting that the 3-D effects are 

not only real in the Monte Carlo but that the mean-field is a good approximation with 

these parameters. In Section [7] we discuss what this expansion really means. 
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Mean Slope Per Segment V.p 




10"' 10" 
Inverse Temperature p) 



Figure 5. This figure shows the mean slope per segment, 5/a, where 5 ~ 0.0098 and a is given 
by Equation 16.21 and that straightness holds for all (3 values. The point at which begins to 
increase with decreasing /3 (Figure[4]) has a mean slope of 35, and, even at the smallest (3 = 10""^, 
the segments have an average angle of 82° with respect to the complex plane. 

6.2. Straightness Holds 
In order to be considered straight enough, we need 



a < ^ = 0.0098. 



(6.3) 



M 1024 

Straightness holds for all f3 values, shown in Figure [5] These conditions hold on average. 
We ignore extreme low-probability cases as not contributing significantly to the statistics. 
These straightness constraints, coupled with filaments having no attractive interactions, 
mean hairpins (kinks or, in quantum terminology, instantons) do not occur. 

Aside from straightness constraints, the reader might question whether allowing vortex 
filaments to entangle violates the model's assumptions. While this is a valid concern, it is 
not an assumption of the model. Since this fluid is almost-everywhere invicid, it contains 
asymptotically small regions of non-zero viscosity, and, while a totally invicid fluid cannot 
allow vortices to change topology (cross over each other) from state to state, an almost- 
everywhere invicid fluid allows vortex reconnections and cross-overs to occur due to 
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microscopic viscous effects. Therefore, in our simulations vortices are allowed to cross 

one-another. We are not claiming to model vortex reconnection, which is a mysterious 

process, but only the before and after effects of it. 

Concerning the question of how vortices can cross one-another and still remain nearly 

parallel, we point to the extremely high-density (tiny value of R^) which allows even the 

straightest filaments to entangle. 

7. Discussion 

We have shown that with increasing (3 there is a "transition" from decreasing R^, the 
2D behavior, to increasing i?^, the 3D behavior. We put this word "transition" in quotes 
because at present we have no proof for or against this being a phase transition. The free 
energy function that we derive is smooth for all positive (3. However, this is not proof 
against there being a phase transition in the original system. At this point it is a subject 
for future research. 

The decreasing-/?, i?-expansion suggests that, by adding degrees of freedom to the 
2D model to make it a quasi-2D model, we add a mechanism for the vortices to resist 
confinement through entropic effects. As (3 decreases past the "transition" , the system's 
i?^ goes from being the result of interaction-versus-angular momentum competition to 
an entropy-versus-angular momentum competition. In the 2D system this 3rd dimension 
entropy is not there. Although there is another kind of entropy in the 2D model that will 
slow the compression as /3 ^ to a constant value, it is not enough to cause an increase 
in the system size. In the 3D system, the degrees of freedom are exponentially greater, 
which causes the expansion seen in Figure (4] 
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8. Conclusion 

Statistical mechanics provides a way to model transitions in decaying turbulence when 
there is a separation of time scales. We conclude from our results that the transfer of the 
model size from an interaction-angular momentum competition to an entropy-angular 
momentum competition is a kind of transition to turbulence purely due to 3D effects. 
Whether it is a phase transition is unknown, but it shows that 3D effects do become 
significant in determining the length scale of the turbulent system, a useful result for 
experiments in tight confinement of rotating nearly-invicid fluids. 

The key role of angular momentum in the derivation of the expression for and 
the subsequent Monte Carlo validation cannot be over-emphasised. At positive /3 and 
a, quasi-2D vortices of the same sign effectively repulse one-another. Therefore, in the 
unbounded plane, they would fly off to infinity without the angular momentum constraint 
unless hemmed in by an infinite expanse of vortices, the case of periodic boundaries, or 
walls. As we mentioned in the introduction, walls and periodic boundaries enforce an 
artificial length scale on the vortices, while angular momentum constraints allow that 
length scale to be found naturally. 

In experiments and simulations of experiments, it is often reasonable to model walls 
present in the experimental setup. However, in oceans and atmospheres and in stars, 
there are no walls, and, at scales where Coriolis effects are small, it does not make 
sense to use a spherical domain unless the simulation is of extremely high resolution. 
The proper regime for applications of unbounded plane, conserved angular momentum 
simulations i s in sm all Rossby number regions where small-scale rotation occurs such as 



Julien et al 



(|1996r ) have modeled in their astrophysical simulations and found arrays of 
like-sign nearly parallel vortex filaments. To this niche we have contributed a previously 
unseen entropic transfer to turbulence in rotating, almost-everywhere ideal fluids. 



Trapped slender vortex filaments 
Acknowledgment s 

This work is supported by ARO grant W911NF-05-1-0001 and DOE grant 

DE-FG02-04ER25616. 



30 T. D. Andersen and C. C. Lim 

REFERENCES 

A. A. Abrikosov. On the magnetic properties of superconductors of the second group. Soviet 

Physics JETP, 5(6):1442-52, 1957. 
V. Berdichevsky. On statistical mechanics of vortex lines. Int. J. Eng. Set., 40:123-129, 2002. 
T. H. Berlin and M. Kac. The spherical model of a ferromagnet. Phys. Rev., 86(6):821, 1952. 
A. J. Callegari and L. Ting. Motion of a curved vortex filament with decaying vortical core 

and axial velocity. SIAM J. on Applied Math, 35(1):148-175, 1978. 
D. M. Ceperley. Path integrals in the theory of condensed helium. Rev. o. Mod. Phys., 67:279, 

April 1995. 

A. J. Chorin. Vorticity and Turbulence. Springer- Verlag, New York, 1994. 

M. T. DiBattista and A. Majda. Equilibrium statistical predictions for baroclinic vortices: The 

role of angular momentum. Th. and Camp. Fluid Dyn., 14:293, 2001. 
R. P. Feynman and J. W. Wheeler. Space-time approach to non-relativistic quantum 

mechanics. Rev. o. Mod. Phys., 20:367, 1948. 
H. Hasimoto. A soliton on a vortex filament. J. Fluid Mech., 51:472, 1972. 
K. Julien, S. Legg, J. McWilliams, and J. Werne. Rapidly rotating turbulent rayleigh-benard 

convections. J. Fluid Mech., 332, 1996. 
R. Klein and A. J. Majda. Self-stretching of a perturbed vortex filament i: The asymptotic 

equation for deviations from a straight line. Physica D, 49:323, 1991. 
R. Klein, A. Majda, and K. Damodaran. Simplified equation for the interaction of nearly 

parallel vortex filaments. J. Fluid Mech., 288:201-48, 1995. 
C. C. Lim. Phase transitions to super-rotation in a coupled barotropic fluid - rotating sphere 

system. Moscow, Aug 2006. lUTAM, Springer- Verlag. 
C. C. Lim and S. M. Assad. Self-containment radius for rotating planar flows, single-signed 

vortex gas and electron plasma. R & C Dynamics, 10:240-54, 2005. 
C. C. Lim and J. Nebus. Vorticity Statistical Mechanics and Monte-Carlo Simulations. 

Springer, New York, 2006. 
P-L. Lions and A. J. Majda. Equilibrium statistical theory for nearly parallel vortex filaments. 

pages 76-142. CPAM, 2000. 



Trapped slender vortex filaments 31 
A. Majda and X. Wang. Non-linear Dynamics and Statistical Theories for Basic Geophysical 

Flows. Cambridge UP, Cambridge, 2006. 
L. Onsager. Statistical hydrodynamics. Nuovo Cimento Suppl., 6:279-87, 1949. 
L. Ting and R. Klein. Viscous Vortical Flows, volume 374 of Lecture Notes in Physics. 

Springer, Berlin, 1991. 
A. Zee. Quantum Field Theory in a Nutshell. Princeton UP, Princeton, 2003. 



Amplitude Per Segment Vs. Segment Height 



Amplitude Per Segment 0(8 ) 
' Segment Height 0(8) 



3 



-3 



10" 



10"' 10^ 
Inverse Temperatuer 



10' 



